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As the number of legs n of an n-leg, t — J ladder increases, density matrix renormalization group 
calculations have shown that the doped state tends to be characterized by a static array of domain 
walls and that pairing correlations are suppressed. Here we present results for a t — t' — J model in 
which a diagonal, single particle, next-near-neighbor hopping t' is introduced. We find that this can 
suppress the formation of stripes and, for t' positive, enhance the d x 2_ y 2-\ike pairing correlations. 
The effect of if > is to cause the stripes to evaporate into pairs and for t' < to evaporate into 
quasi-particles. Results for n — 4 and 6-leg ladders are discussed. 



PACS Numbers: 74.20. Mn, 71.10.Fd, 71.10.Pm 

Neutron scattering experiments on Lai.6- x Ndo.4 
Sr^CuO^ show evidence of a competition between static 
(quasi-static) stripes and superconductivity [fij. Here 
the stripes consist of (1,0) domain walls of holes sepa- 
rating 7r-phase shifted, antiferromagnetic regions. For 
x = 0.12 (x w 1/8), the intensity of the charge and 
spin superlattice peaks is largest and T c is less than 5K. 
As x deviates from this value, the relative intensity of 
the magnetic superlattice peaks decrease and the super- 
conducting transition temperature T c increases. High 
field magnetization studies H provide evidence that in 
this material superconducting can coexist with static (or 
quasi-static) stripe order. However, the fact that T c is a 
minimum where the superlattice peaks are most intense 
suggests that static stripe order competes with supercon- 
ductivity. 

We are interested in understanding whether at — J-like 
model can exhibit this type of behavior. In studies of n- 
leg, t — J ladders we have previously found evidence for 
stripe formation. In particular, for n = 3 and 4 legs we 
have found evidence for both stripes and pairing H|J. 
These systems have open boundary conditions. How- 
ever, in wider ladders (n = 6 and n — 8) with cylindrical 
boundary conditions, where the stripes close on them- 
selves rather than having free ends, the stripes appeared 
to be more static and the pairing correlations were found 
to be suppressed . This suppression of the pairing cor- 
relations was also observed when an external potential 
was applied to further pin the stripes. 

If the formation of static stripes could be suppressed, 
one might hope to find enhanced pairing correlations. It 
is not clear whether the complete elimination of stripes or 
only a slight dcstabilization would be more favorable to 
pairing correlations. We have been investigating various 
interaction terms which could destabilize stripes. Here 
we focus on the effect of a next-near-neighbor diago- 
nal hopping t' . Effective hopping parameters have been 
evaluated from band structure calculations and finite 
CuO cluster calculations. For the hole-doped cuprates 
t' is found to be negative while for the electron-doped 



cuprates it is positive. Both t' and the one-electron hop- 
ping t", which connects next-near-neighbor sites along 
the (0,1) or (1,0) axis, have been used in t — t' — t" — J 
models to fit ARPES data Q. In addition, Lanczos cal- 
culations by Tohyama and Maekawa [pj on t — t' — J clus- 
ters and Monte Carlo calculations on t — t' Hubbard 
lattices show that t' > tends to stabilize the commensu- 
rate (tt, 7t) antiferromagnetic correlations. Recently, ex- 
act diagonalization and density-matrix renormalization 
group (DMRG) calculations on small clusters and four- 
leg ladders have found that t' < destabilizes stripes |j| . 
Furthermore, it was concluded that a small positive t' 
did not destabilize the stripes on these systems. 

Here we will consider the effect of t' on both open four 
leg and cylindrical six leg ladders. In addition to consid- 
ering the affect of t' on stripe stability, we will measure 
its affect on pairing correlations. We find that stripes 
are destabilized for either sign of t', and that pairing is 
suppressed for t' < 0, and enhanced for t' > 0. This lat- 
ter effect is surprising, since superconducting transition 
temperatures are generally higher for hole doped cuprates 
(f < 0) than for electron doped (*' > 0) . 

The t — t' — J Hamiltonian which we will study has the 
form 

+ J Y,(ij) 0i ■ Sj - \ nrij) (1) 

Here (ij) are near-neighbor sites, (ij)' are diagonal next- 
near-neighbor sites, Si = \c~l s cr ss >Ci S , m = c^c^ + c^c,j., 
and cf s (ci S ) creates (destroys) an electron of spin s at site 
i. No double occupancy is allowed. We will use DMRG 
calculations to explore the charge, spin, and pairing cor- 
relations on doped four and six leg ladders. The calcu- 
lations reported below keep up to 1200 states per block, 
with truncation errors of about 10 -4 , and from six to ten 
finite system sweeps. We have checked the inclusion of t' 
in our program by comparing the results for the rung hole 
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density on a 14 x 4 system with the results of Tohyama, 
et. al. ||; precise agreement was found. 

In a previous study of the 4-leg t — J ladder, we found 
that four-hole diagonal domain walls formed as the dop- 
ing increased. In Figure 1(a) and (b) we show the rung 
density 



For this same 12 x 4 lattice, we have studied the pair- 
field correlation function 



(2) 



i=l 



versus I for J/t = 0.35 on a 12 x 4 lattice with 8 holes 
and open boundary conditions. For t' = 0, we clearly see 
the formation of two domain walls, signaled by two broad 
peaks in (n r (£)). As t'/t is varied, one clearly sees that 
the static domain wall structure is suppressed for either 
sign of t'. 
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FIG. 1. (a) Hole density per rung for a 12 x 4 system 
with 8 holes, J/t = 0.35 and open boundary conditions, with 
i <= 0. (b) Same as in (a), but with t' >= 0. (c) and (d) 
d-wave pairing correlations for the systems shown in (a) and 
(b), respectively. 



D{t) = (A i+i Af) 



(3) 



with Af a pair creation operator which creates a singlet 
d x 2_ y 2 pair centered on the i th site of the second leg. 
Figure 1(b) shows a plot of D(£) versus t for the 4 x 12 
ladder for J/t = 0.35 with 8 holes and various values of 
t'/t. As t'/t initially increases, the pairing correlations 
are enhanced but as t'/t becomes greater than ~ 0.3, 
they are suppressed. They are suppressed for t' negative, 
with very strong suppression occuring for t' < —0.2. 

Results for the charge density and spin structure of 
a 12 x 6 lattice with J/t = 0.5 and 8 holes are shown 
in Figure 2. Here we have taken cylindrical boundary 
conditions, i.e. periodic in the y-direction, open in the 
x-direction. In this case, for t'/t — 0, the holes form two 
transverse domains each containing 4 holes. The 7r-phase 
shifted antiferromagnetic regions which are separated by 
these domains are clearly visible in Figure 2 for t' = 0. 
The DMRG calculation has selected a particular spin or- 
der, breaking symmetry; as the number of states kept 
per block increases, the magnitude of this spin order de- 
creases, and the exact ground state would have no net 
spin on any site. However, here the spin order serves 
to illustrate the underlying spin correlations in the exact 
ground state, which we expect to be a superposition of 
the broken symmetry state rotated to all possible direc- 
tions. 
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FIG. 2. Hole and spin densities on 12 x 6 systems with 
cylindrical boundary conditions. The hole density is propor- 
tional to the diameter of the circles, according to the indicated 
scale, and similarly the length of the arrows gives the expec- 
tation value of S z . 
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As t'/t increases, we again see a suppression of the 
charge order and in addition the 7r-phase shifted anti- 
ferromagnetic regions disappear. This is also true for t' 
negative. For t' = 0.3, we see that Neel spin order, with- 
out any 7r phase shifts, is now the broken symmetry state. 
As previously noted, Lanczos and Monte Carlo B cal- 
culations indicated that a positive t' tended to stabilize 
the commensurate (ir, tt) antiferromagnetic correlations, 
which is consistent with our results. 

The rung density shown in Figure 3(a) provides a more 
quantitative display of the suppression of the charge do- 
mains walls. In this finite magnitude of t' seems 
to be necessary to substantially reduce the charge density 
structure. The domain walls in L x 6 ladders at t' = are 
stable bound states of two hole pairs, and a finite change 
in the parameters of the systems is needed to break them 
up. We believe that L x 6 cylindrical systems have un- 
usually stable domain walls, and that more generally a 
smaller value of \t'\ would destabilize the stripes. Here, 
we see that the stripes are suppressed for t' = 0.2, and 
completely destabilized for t' = 0.3. 

Figures 3(c) and (d) show the pair-field correlations 
D(l) versus £ for various values oit'/t for the 12 x 6 lad- 
der. We see when the stripes are weakened by a positive 
t' , pairing correlations are strongly enhanced. The opti- 
mal t' appears to be near t' = 0.2. Pairing is once again 
suppressed for negative t', even when the domain walls 
are destabilized. 

From a weak coupling point of view, our results on the 
effect of t' on pairing are surprising. In weak coupling, 
the effect of t' < is to shift the van Hove singularity in 
the density of states away from half-filling, so that the 
singularity may occur near the Fermi level in a doped 
system. Thus, one might have expected to find an en- 
hancement in pairing for t' < 0. However, in the t — J 
model, we find a suppression of the pairing. In strong 
coupling, one can understand this effect. Consider a pair 
of holes, and imagine we fix one hole and let the other 
hole hop around it. Consider the phase of the wavefunc- 
tion of the second hole on the four sites next to the first 
hole. It appears that t' < will directly favor a +-+- 
d-wetve phase pattern as the second hole hops around 
the first, whereas t' > would favor the s-wave pat- 
tern ++++ [flOfl . However, the actual phase of a pair 
is a relative phase between a system with N holes and 
one with N + 2 holes. If one considers a 2 x 2 t- J sys- 
tem, one finds that the 2-hole ground state has s-wave 
rotational symmetry, whereas the undoped state has d- 
wave rotational symmetry [|TT|-|l3| . The d-wave nature of 
the pairing comes from the difference in these rotational 
symmetries. Consequently, t' < 0, by suppressing the 2- 
hole H — H — h pattern, actually suppresses d-wave pairing, 
while t' > can enhance it. 
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FIG. 3. (a) Density of holes per rung for the 12 x 6 ladder 
systems shown in Fig. 2. (b) d-wave pairing correlations for 
the same systems. 

Consider the 2x2 system The energy of the 

undoped system is independent of t'; we find £"(0) = 
—3 J. The energy of the one hole system depends only 
weakly on t'\ for t' small, we find £7(1) = - J - 1/2(J 2 + 
12< 2 + 4Ji' + 4i' 2 ) 1 / 2 . For J = 0.35, t = 1, this varies with 
t' as £7 w -2.09087 - Q.1005t'. The energy of the two- 
hole system, in contrast, depends strongly on t': £7(2) = 
- J/2 - 1' - (32t 2 + ( J + 2f ) 2 ) 1 ' 2 . The pair binding energy 
is defined as 



E b = 2£7(1) - £7(2) - £7(0). 



(4) 



The dependence of the pair binding energy on t' is domi- 
nated by £7(2), and we find that t' > Q strongly enhances 
the pair binding. 

On larger systems, the detailed energetics are more 
complex, but a similar effect occurs. In Fig. 4, we show 
the energy per hole of several systems as a function of 
t' Q. The systems allow us to compare the stability 
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of paired states, striped states, and states with isolated 
holes. The first system is a single hole in an 8 x 8 open sys- 
tem, with a staggered antiferromagnetic field of strength 
0.1 on the edges to approximate the magnetic coupling 
to the rest of the system, which is assumed to be un- 
doped. The second system is similar, but has two holes. 
We plot the energy difference between these systems and 
the same system without holes, divided by the number 
of holes [jnj. The third system is a 16 x 6 system, with 
open boundary conditions, and staggered fields of mag- 
nitude 0.1 with a 7r phase shift applied on the first and 
last chain. These boundary conditions favor the devel- 
opment of a stripe down the center of the ladder. Then 
we subtract the energy of an undoped 16 x 6 system, also 
with staggered fields, but without the phase shift Q. 
We expect that finite size effects are not neglible, and 
these could shift the striped phase curve relative to the 
other two curves. However, we believe the general trends 
are reliable. That is, the striped system is lowest in en- 
ergy near t' = 0, but becomes unstable as t' becomes 
less than —0.1, or as t' increases above a value slightly 
greater than 0.0. Thus, the striped region is quite nar- 
row as a function of t' . This conclusion differs somewhat 
from that of Ref. || , where it was found that stripes were 
enhanced for < t' < 0.2, but were suppressed for larger 
values of t'. For positive t', the new stable state has pairs 
of holes, as Tohyama, et. al. f| also found for t' ~ 0.5. 
For t' < —0.1, the near degeneracy between one and two 
holes indicates that the holes are not bound into pairs: 
instead, the stripes break up into quasiparticles. These 
observations are consistent with enhanced pairing corre- 
lations for t' > 0, and suppressed pairing correlations for 
t' < 0. Note that 
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FIG. 4. Energy per hole of various hole configurations, as 
discussed in the text. 

iV-leg, t — J ladders provide perhaps the simplest mod- 
els which exhibit many phenomenologically similar char- 
acteristics to those observed in the cuprates. Here, for 
two different t — t' — J ladders, we find that a diagonal, 
next-near-neighbor hopping suppresses the formation of 



static stripes and that for t' > this can lead to an en- 
hancement of the d x 2_ y 2 pairing correlations, while t' < 
we find suppression of pairing. 
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